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Abstract 

The quark potential model is extended to include the sea quark excitation using the random 
phase approximation (RPA). The effective quark interaction preserves the important Quantum 
Chromodynamics (QCD) properties — chiral symmetry and confinement simultaneously. A pri- 
mary qualitive analysis shows that the vr meson as a well-known typical Goldstone boson and the 
other mesons made up of valence qq quark pair such as the p meson can also be described in this 
extended quark potential model. 
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I. INTRODUCTION 



The Quantum Chromodynamics (QCD) is widely believed to be the correct theory of 
strong interaction nowdays. Due to the complication of non-Abelian gauge theory, in the 
studies of hadron structure, models are often adopted to approximate the strong interaction. 
Two kinds of models are often used, each incorporates several important QCD features. 
The quark potential model based on the confinment quark picture has been impressively 
successful in hadron spectra [^, except for the vr meson. Recently, it was shown that for 
the TT meson, which is a well-known Goldstone boson, the explanation of quark potential 
model, i.e. mesons merely made up of qq pairs of valence quarks, is not well satisfied and the 
important sea quark contribution must be considered 0. On the other hand, the Nambu- 
Jona-Lasinio (NJL) model and several of its extensions 0, ^ can well describe the vr 
meson as a Goldstone boson which is a conclusion of QCD chiral symmetry. The sea quark 
contributions to the meson structure are included in the NJL model by using the Bethe- 
Salpter equation or its equivalence — Random Phase Approximation (RPA). Also ground 
states of other mesons such as p which are made up of simple valence qq quark pair can be 
described by the extensions to the original NJL model. However, the NJL model lacks the 
important QCD property of quark confinement, so the excitation states of mesons cannot 
be accounted in the model, such as tt', p', p", etc. 

In this paper, we try to incorporate the important QCD chiral symmetry into the quark 
potential model by comprising the sea quark contribution. To incorporate the chiral symme- 
try, we start from an effective quark hamiltonian taken from extended NJL models which is 
chiral invariant. To consider the sea quark contribution, we adopt the RPA approximation. 
This derives a hamiltonian with two channels, i.e. the ordinary valence quark channel and a 
new sea quark channel, which is an extension to the quark potential model. This extended 
quark potential model with RPA approximation demands a new potential which couples 
the valence quark channel with the sea quark channel we call coupling potential. Then a 
preliminary numerical analysis is made to illustrate that the Goldstone vr meson and also p 
meson can be well described together in the extended quark potential model. 

In the next section. We obtain the extended quark potential model — the RPA approxi- 
mation is used to derive the two channel hamiltonian and the so called coupling potential. 
In preliminary analysis is made to the vr and p mesons. Finally, some brief summary 
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and discussion are given. 



II. EXTENDED QUARK POTENTIAL MODEL WITH RPA 

In this section, we will incorporate the QCD chiral symmetry into the original quark 
potential model. To avoid the complication of QCD, our start point will be an effective 
quark interaction. A good choice of the effective quark interaction with chiral symmetry is 
of the NJL type. However, despite of much success of the original NJL model, its contact 
quark interaction is too simplified as an effective quark interaction which does not possess 
the important QCD feature of color confinement. The original NJL model has been modified 
by several authors to include the nonlocal effective quark interaction 

Consider a quark nonlocal action where the QCD strong interaction is approximated 
by the quark effective interaction after the complicated degrees of freedom of gluon are 
eliminated as in the extended NJL models. Its general form is 



where K^{x) are the interaction kernel functions and are tensor products of Dirac, flavor 
(in this paper, we will restrict to iso-spin SU{2) symmetry only) and color {SU{3)) matrices 
which represent the quark interaction channels. The quark interaction in this quark action 
will be chosen to have chiral symmetry and with quark confinement in the next section. 

As the action is nonlocal, the exact hamiltonian will be very complicated and also time 
dependent. In the present paper we shall only consider the non-relativistic limit. We can 
adopt a simple, time independent hamiltonian with an instantaneous approximation. After 
approximating the kernel function as 



A 




(1) 



(2) 



the model hamiltonian becomes 




H 



4, 



(3) 
(4) 



(5) 
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where "^{x) = \E'(a3,0) which can be expressed in quark creation and annhilation operators 

v[/(a.) = 5: y [b{k, s)u{k, 8)6^"- + d\k, s)v{k, 8)6-^"-] . (6) 

The hamiltonian from the instantaneous approximation breaks the important Lorentz in- 
variance. We shall only use it to study the hadron structure in rest frame. For a moving 
hadron, however, we shall use the Lorentz boost technique to restore the Lorentz invariance 
in future study. 

Apart from a constant, the hamiltonian (|]) can be written in normal order form as 

H = Ho+:H,: (7) 

where Ho represents the single particle energy of constituent quark and it can be re-written 
formly as 

Ho = J d^x^^ix) [-iA{-V^)a. ■ V + /^^(-V^)] ^{x). (8) 

The function forms of A and B operators should be obtained from the Mean Field Approxi- 
mation (MFA) and their momentum dependence are due to the Lorentz invariance breaking 
from instantaneous approximation. Here for simplification we shall take the A and B op- 
erators as constants. This only introduce one parameter of the constituent quark mass rric- 
The constituent quark single energy is simplified as 

Ho = J d^x^^{x) [-la. ■ V + /3m,] ^{x). (9) 

It is generally accepted that hadrons are excited states of the QCD vacuum. To account 
the sea quark contribution, we consider a meson as a superposition of creation and annihila- 
tion of q-q pairs on the vacuum. Thus the meson states are treated as an excitation modes 
of the vacuum of the RPA type: 

\Q) = |0) , (10) 

where is the excitation operator. For a meson in the rest frame, can be expressed in 
quark and antiquark creation and annihilation operators 

r d^k 1 r 
= J2j (27r)3 2fc0 ^"^^^^ S2)b\ksi)d^{-ks2) 

+y{k, si, S2)d{-ksi)b{ks2)] , (11) 
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where x and y are the x-amphtude and y-amphtude of the RPA wave function respectively 
and the time reversal operators are 

h{ks) = {-iy^^-'b{k-s), (12) 
d{ks) = {-iy/^-'d{k-s). (13) 

The X- and y-amplitudes can be solved from the RPA equation of motion 

(0| [5Q, [H, Qt]] |o) = E (0| [5Q, Q^] |0) . (14) 

After inserting the hamiltonian eq. (^ into the RPA equation of motion, we obtain the 
equations for the x- and y- amplitudes 

[e -2k^)x{k,si,S2) 

{u{ksi)r^v{—ks2)v{—k's4)T^u{k's3)x{k', S3, S4) 
-u{ksi)T,v{-ks2)'ii{k's4)T,v{-k's3)y{k', S3, S4)} 

{u{ksi)Tf:u{k's3)v{—k's4)Ti:v{—ks2)x{k', S3, S4) 
-u{ksi)T^v{-k's3)'a{k's4)T^v{-ks2)y{k', S3, S4)} , (15) 

(E + 2A;°)y(fc,si,S2) 



and 



^v{~ksi)T ^u{ks2)v{—k' S4)T ^u{k' S3)x{k' , S3, S4) 
-lj{-ksi)T,u{ks2)'a{k's4)T,v{-k's3)y{k', S3, S4)} 



S3S4 



^v{—ksi)T ^u{k' S3)v{—k' S4)T eu{ks2)x{k' , S3, S4) 

-v{-ksi)Tev{-k' S3)u{k' S4)T^u{ks2)y{k' , S3, S4)} , (16) 

where u and v are the time reversal of spinors u and v respectively. In each equation, the 
first summation (direct term) can be absorbed into the second summation (exchange term) 
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by Fierz rearrangement, then equations briefly become 



and 



E - 2rj x(fc,si,S2) 

{u{ksi)V ^u{k' S3)v{—k' S4)V ^v{—ks2)x{k' , S3, S4) 
-u{ksi)T^v{-k's-i)U{k's4)T^v{-ks2)y{k\ S3, S4)| , (17) 

(E + 2A;°)y(fc,si,S2) 

— fcSi)reM(fc's3)t;( — fc's4)reU(fcS2)a;(fc', S3, S4) 

-^-ks{)V,v{-k's^)u{k's^)V,u{ks2)y{k', S3, S4)} . (18) 
Now We introduce two effective potentials U and V, such that the matrix elements are 
(fc, si, S2 I f/ I k\ S3, S4) = ^^(^ - 

M(fcSi)re'u(fc's3)t;( — fc's4)rel'( — fcS2), (19) 

(fc,si,S2 I V I A;',S3,S4) = 

M(fcsi)r,^y(-fc's3)S(fc's4)r,?;(-fcs2). (20) 

From the basic sjTiimetries of parity, charge conjugation and time reversal of quark in- 
teraction, and the definition eqs. ([I9|) and (pOD , one can obtain the hermicity of effective 



potentials U and V 

= U, (21) 
= V, (22) 



and also 



(fc, Si,S2 I U I fc', S3,S4)* 

(_l)l/2-.5i(-_]^^l/2-S2^_-|^^l/2-S3^_]^^l/2-S4 

(fc, -S2, -Si I t/ I fc', -S4, -S3). (23) 



Finally, the RPA equations are written as 

- 2k°) x{k,si,S2) 



J2 J 7K-^ [(^' ■^i' ^2\U \ k', S3, S4)x{k', S3, S4) 



S3S4 



(27r) 



and 



-{k, si, S2\V \ k\ S3, Si)y{k', S3, s^)] , (24) 

(E + 2A;0)y(fc,si,S2) 
= / TTTII [(^' ■^i' ■52 I ^ I k', S3, Si)x{k', S3, Si) 

S3S4 v^^/ 

-(A;, si, S2\U \ k', S3, Si)y{k', S3, S4)] . (25) 
In the above RPA equations, the potential U interacts between the same RPA amphtudes 



(x-amphtude in eq. (|2J) and y-amphtude in eq. (^Sf )). The x-amphtude represents the con- 
tribution of valence quarks, that is just the wave-function of valence quarks in the potential 
model. We conclude that the U potential is just the quark potential in the potential model. 
On the other hand, potential V interacts between the different RPA amplitudes. This is a 
new potential and we call it the coupling potential. 

The above RPA approximation may be considered as an extension to the potential model. 
The extended potential model includes two channels — the ordinary valence quark chan- 
nel and an extended sea quark channel. The valence quark channel which represents the 
contribution of valence quarks is described by the x-amplitudes. The sea quark channel 
representing the contribution of sea quarks is described by the y-amplitudes. The ordinary 
quark potential U interacts within both individual channel. The new coupling potential V 
couples the two channels together. 

The potentials U and V will be evaluated from the kernels in action (|lD. Following ref. 0, 
from the spinor structure of Lorentz invariant, the kernels can be divided into five types, 
i.e., scalar {S), pseudo-scalar (P), vector {V), axial- vector {A), and tensor (T). We write 



(26) 

In the non-relativistic limit k, k' 0, only the following terms survive in the potentials 

U = -S + V + Aa^ - (T2 + T(Ti- (72, (27) 
V = P - VcTi -(72- A- Ta-i ■ 0-2. (28) 

where cr's are the pauli spin matrices. 



III. PRIMARY APPLICATION 



In addition to the Lorentz invariance of quark interaction kernels, in the hmit of vanishing 
current quark masses, the kernels must also have the chiral symmetry of QCD. The general 
quark interaction which is compatible with the chiral symmetry can be found in ref. 0. 
Here we simply choose the form of effective quark interaction H4 as a vector interaction like 
the one gluon exchange interaction 

H4 = \j d^xd'yK{x-y)^{x)^^X'^^i^{x)^{y)^^\'^^{y). (29) 

Following the quark potential model |^ , the effective quark interaction kernel K may be 
chosen as a confinement interaction plus the one gluon exchange interaction 

Kir) = -^^br + \^. (30) 

Our choice of confinement interaction differs from ref. in that ours is a vector 7^^ ® 7^ 
interaction. This is because from the effective interaction kernel eq. (pB]), according to ref. [0, 
the effective qq interaction is 

U = S + V - Acri- (72 + Ta-i- (T2. (31) 

Notice that scalar interaction 1 1 changes sign from the qq (eq. (0)) to qq (eq. (|3T|)). This 
seems to us that one can not get the confinement for both mesons and baryons simutaneously 
from a unified interaction kernel. 

To actually perform calculations with these potentials, we would like to diagonalize the 
hamiltonian matrix in the harmonic-oscillator basis following ref. However, in this 
primary analysis, we simply approximate the meson ground states using the ground state 
wave function of harmonic-oscillator 

x{k,si,S2) = X^ooo(A;)(l/2sil/2s2 | SMg), (32) 
yik,s,,S2) = y^ooo(fc)(l/2sil/2s2 | SMs), (33) 

where SMs is the spin of the meson, and 

^ooo(A:) = y^-^'^'^"^- (34) 
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From the above simphfication and the effective hamiltonian eq. (]29|), the RPA equations 
reduce to 



(35) 



A 


B 




X 


= m 


X 


-B 


-A 




Y 




Y 



where 



A = 2{k') + {K), 



(36) 
(37) 



and 



(cTi - era) = 2^(^ + 1) -3. 

In above formula, (...) takes the average value. In nonrelativistic limit k 
= ^k'^ + ml ~ rric, so we have approximately 

A ^ 2mc + {K). 

By introducing 

a = X + Y, 

/? = X - y, 



(38) 

0, hence 

(39) 



(40) 
(41) 



the RPA equations become 



{A + B)a = mf3, 
{A-B)p = ma. 



(42) 
(43) 



For a preliminary analysis, we will only consider two typical mesons — the well-known 
Goldstone boson tt and the vector meson p which is well described in the original quark 
potential model as qq valence quark pair. For vr, its orbital angular momentum L = 0, spin 

= 0, we have 



2{mc + 2{K))a^ = m^P^, 
2{mc - {K))p^ = m^a^. 



(44) 
(45) 
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The chiral symmetry gives a relation between constituent quark mass rric and average 
value of effective quark interaction (K) . From MFA 



rrir 



rUg + E, (46) 



where rriq is the current quark mass, S is the Hatree-Fock quark self-energy. In the chiral 

limit rriq — > 0, the pion will be an exact massless Goldstone boson, i.e., 0. From the 
equations of pion, eqs. ( PD and p5|), we obtain in the chiral limit where rriq = 0, = 0, 

E= (K). (47) 

Finally 

rric = mq + (K). (48) 
Back to equations of vr meson, we obtain 



m„ = 2^'mq{3mc-2mq), (49) 

(50) 



"tt _ mq 



Pt^ Y Snic — 2mq 
For p meson, its L = and 5' = 1, then 

2mcap = ruppp, (51) 
2{mc + {K))f3p = nipap. (52) 

We obtain 

rup = 2^mc(2mc - m^), (53) 
Op ^ \2m^-mq_ 

(3p V 

In the table |, we give two sets of numerical calculations. In the first calculation (A), the 
quark masses {rriq and rric) are taken to be general accepted values. In the second calculation 
(B), the quark masses are chosen to better fit the experimental data of meson masses (m^ 
and rrip). We believe that both sets of results are reasonable as a primary calculation. 
The amplitude ratio Yj^^p/X^^^p determines the sea quark contribution vs. the valence quark 
contribution 0. Since ~ the sea quark contribution is essential to vr meson. For p 
meson, Xp ^ Yp, the sea quark contribution is small. 
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TABLE I: Numerical calculations (all mass units are MeV). 





ruq 


rUc 




nip 






A 


5.5 


300 


110 


845 


-1.2 


5.9 


B 


6 


274 


139 


771 


-1.2 


5.9 



Since m^r S> m^, we have approximately 

= 2-\/3y^m^m^, (55) 

mp = 2v^mc, (57) 



Then we find a mass relation 



ml = VTSmqUip (59) 



The experimental data m.,^ — 139MeV, rup — 770 MeV combining with the general accepted 
estimation rUq — 5.5MeV give 

,2 



= 4.6, (60) 

niqmp 



while ^/TS = 4.2. 



IV. SUMMARY AND DISCUSSION 



To take into account of the sea quark contribution to mesons more properly, the original 
quark potential model is extended which now includes a new sea quark channel similiar to 
the RPA approximation. A new coupling potential is needed to couple the original valence 
quark channel and the sea quark channel. Through an effective quark interaction which is 
also chiral invariant as demanded by QCD, the new potential is related to the conventional 
quark potential in original quark potential model. The constituent quark mass is related to 
the quark potential also from the chiral symmetry. The primary calculation connected the 
masses of mesons to masses of quarks (constituent and valence) . The sea quark contribution 
is important to the tt meson which is a well-known Goldstone boson. Other mesons, such as 
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vector p meson, where the valence quark contribution is dominant, can also be well described 
by the same extended quark potential model. 

Finally, we shall indicate that in this primary analysis, under the non-relativistic approx- 
imation, many important features of quark interaction such as the spin-orbit interaction, 
the spin-spin interaction and tesor interaction are neglected. Those must be considered in 
the complete calculation of meson spectroscopy in the future. 
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